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ABSTRACT 



We report on spectroscopic and photometric observations through transits of the exoplanets WASP-22b and WASP-26b, intended to 
determine the systems' spin-orbit angles. We combine these data with existing data to refine the system parameters. We measure a 
sky-projected spin-orbit angle of 22+16° for WASP-22b, showing the planet's orbit to be prograde and, perhaps, slightly misaligned. 
We do not detect the Rossiter-McLaughlin effect of WASP-26b due to its low amplitude and observation noise. We place 3-<r upper 
limits on orbital eccentricity of 0.063 for WASP-22b and 0.050 for WASP- 26b. After refining the drift in the systemic velocity of 
WASP-22 found by Maxted et al. (2010, AJ, 140, 2007), we find the third body in the system to have a minimum-mass of 5.3 ± 0.3 
Mj U p (a 3 /5 AU) 2 , where a 3 is the orbital distance of the third body. 

Key words, binaries: eclipsing - planetary systems - stars: individual: WASP-22, WASP-26 



1. Introduction 

By taking spectra of a star whilst an exoplanet transits across it 
we can measure the sky-projected obliquity A, where obliquity 
is the angle between the stellar rotation axis and the planetary 
orbital axis. As the planet obscures a portion of the rotating star 
it causes a distortion of the observed stellar line profile, which 
manifests as an anomalous radial-velocity ( RV) signature known 
as the Rossiter-McLaugh lin (RM) effect dHold Il893t iRossited 
11924 lMcLaughlinlll924|) . The shape of the RM effect is sen- 
sitive to the path a planet takes across the disc of a star relative 
to the stellar spin axis, thus by observing it we can determine A. 
If we can determine the stellar inclination with respect to the sky 
plane / then we can determine the system's true obliquity X V. For 
example, the light curve of a photospherically active star may 
exhibit rotational modulation, from which the true stellar rota- 
tion velocity can be determined. This can be compared with the 
spectroscopically-determined, sk y-projected rotation velocity to 
get th e stellar inclination (e.g. IWinn et all l2007t iHellier etldl 
1201 11) . 

It is thought that the obliquity of a short-period, giant planet 
is indicative of the manner by which it arrived in its current or- 



* Based on observations made with the HARPS spectrograph on 
the 3.6-m ESO telescope (proposal 085.C-0393), the 0.6-m Belgian 
TRAPPIST telescope, and the CORALIE spectrograph and the Euler 
camera on the 1.2-m Euler Swiss telescope, all at the ESO La Silla 
Observatory, Chile. 

** The photometric time-series and radial velocity data used in 
this work are only available in electronic form at the CDS 
via anonymous ftp to cdsarc.u-strasbg.fr (130.79.128.5) or via 
http://cdsweb.u-strasbg.fr/cgi-bin/qcat?J/A+A/ 



bit from farther out, where it presumably formed. As the angular 
momenta of a star and its planet-forming disc both derive from 
that of their parent molecular cloud, stellar spin and planetary or- 
bital axes are expected to at least initially be aligned. Migration 
via tidal interaction with the gas disc is expected to preserve this 
initia l spin-orbit alignment dLin et al. l fl996r. lM arzari & Nelson! 
2009), but a significant fraction of the few dozen systems so 
far measured!] are misal igned and some are even ret rograde (e.g. 
seejHebrard et al. 201 1, and the references therein). Triau d et alj 
(1201 Oh interpreted this as indicating that some or all close- 
in planets arrive in their orbits by a mixture of: planet-planet 
scattering and the Kozai mechanism, which can drive planets 
into eccentric, misaligned orbits; and tidal friction, which cir- 
cularises and shortens orbits (e.g. , Fabrvckv & Tremaine 2007; 
i Naeasawa et al.l2008tlMatsumura et al.ll2010l:lNaoz et al.l201 ll) 
iMorton & JohnsonT i 201 lb found the current spin-orbit distribu- 
tion to be well described by migration via a combination of such 
alignment-preserving and misaligning mechanisms. 

IWinn et al.l (1201 Oal) and lSchlaufmanl (2010) showed that the 
stars in misaligned systems tend to be hot (T e s > 6250 K). 
Win n et al.l (1201 Oal) suggested that this indicates either that plan- 
etary formation and migration mechanisms depend on stellar 
mass, or that cooler stars, with their larger convective zones, are 
able to more effectively realign orbits via tidal dissipation. The 
number of possible exceptions to this trend is increasing, such 
as HAT-P-9 with T eS = 6350 K and A - -16 + 8°, though rather 
than having been reali gned, such syst ems may simply have never 
been very misaligned (Moutou et al. 201 1). 



1 Rene Heller maintains a list of measurements and references at 
http: //www. aip.de/People/rheller/content/main_spinorbit .html 
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In addition to the scattering route, it has been suggested that 
misaligned systems can ar ise when stars are misaligned with 
their discs (lLai et al.ll2.01 ll) or when discs are misaligned with 
their stars (Bat e etalJl2010h . 

In this paper we present obser vations of spectro scopic and 
photometric transits of WASP-22b dMaxted et alJ 2010. hereafter 
M10) and WASP-26b dSmallev et al.ll2010l hereafter S10). from 
which we attempt to determine the obliquities and refine the pa- 
rameters of the systems. 

WASP-22b is a 0.59-M Jllp planet that transits its solar metal- 
licity ([Fe/H] = 0.05 ± 0.08) GOV, V=12.0 host star every 3.53 
days. M10 found a linear trend in the RV measurements of 
WASP-22 suggestive of either a second planet, a low-mass M- 
dwarf or a white dwarf companion. 

WASP-26b is a 1.03-Mj up planet that transits its solar metal- 
licity ([Fe/H] = -0.02 + 0.09) GOV, V=11.3 host star every 2.76 
days. S10 found WASP-26 to be a visual double with a K-type 
star, separated in the sky-plane by 3 800 AU, and used it to infer 
a common age of 6 ± 2 Gyr. 



2. Observations 

2. 1 . Spectroscopic transits and orbits 

A spectroscopic transit comprises a large number of RVs taken 
in quick succession, during which the target star may have a spe- 
cific activity level that could bias the measured systemic velocity 
and hence the determination of A. Typically, the stellar rotation 
period is much longer than the orbital period. So, to first order, 
stellar activity would manifest as a slow variation in the appar- 
ent RV of the system's centre of mass, which is negligibly small 
on the time-scale of a transit. To mitigate against the potential 
effects of activity, we obtained spectra of the two host stars out- 
side of transit on the nights of the measured RM effects, and to 
refine the planets' masses and orbital eccentricities we obtained 
spectra across the full orbits. 

We obtained spectra of WASP-22 and WASP-26 using the 
HARPS spectrograph on the ESO 3.6-m telescope at La Silla. 
Calibration spectra of the thorium-argon lamp were taken at the 
start of each night, thus avoiding contamination of the stellar 
spectra by the lamp. This is made possible by the stability of 
HARPS, which is at the 1 m s _1 night 1 level. 

On 2009 Nov 19, we obtained 20 spectra of WASP-22 
through a transit with HARPS. The exposures were 15-20 min 
and the signal-to-noise ratio (SNR) per pixel at 550 nm was in 
the range 16-34, with 20 being typical. During the sequence, 
the airmass of the target decreased from 1.11 to 1.01 and then 
increased to 1.70. We sampled the full 3.53-day orbit of WASP- 
22b by obtaining a pair of spectra per night, well-spaced in time, 
on the night of the transit and on each of the preceeding three 
nights. With exposures of 30 min for these eight spectra, the 
SNR was 24-48. 

On 2010 Sep 12, we obtained 24 spectra of WASP-26 
through a transit with HARPS. The exposures were 10-15 min 
and the SNR per pixel at 550 nm was in the range 10-32, with 
22 being typical. During the sequence, the airmass of the target 
decreased from 1.62 to 1.03 and then increased to 1.08. As the 
focus was improved at BJD (UTC) = 2 455 45 1 .641 the SNR and 
measurement precision improved, despite a decrease in the ex- 
posure time from 15 to 10 min. We sampled the full 2.76-day 
orbit of WASP-26b by obtaining six spectra over a period of five 
days around the transit. With exposures of 15-20 min, the SNR 
of these spectra was 26-40. 



After removing the instrumental blaze function dTriaud et al.l 

1201 Oh . radial-velocity (RV) measurements were computed for 
each star by wei ghted cross-correlation with a numeri cal G2- 
spectral template dBaranne et al.ll 1996k iPepe et al.ll2005h . These 
are available online only. We ascribe the fact that the first five 
RVs in the transit sequence of WASP-26 are lower precision than 
the rest to higher airmass and poorer focus. 

We also incorporated in our analysis the 37 CORALIE RVs 
and 6 HARPS spectra of WASP-22 reported in Ml and the 16 
CORALIE RVs of WASP-26 reported in S10. For consistency 
with the new measurements, we recalculated RVs from the ex- 
isting HARPS spectra of WASP-22 after removal of the instru- 
mental blaze function, which had not been done previously. We 
present these recomputed RVs along with the new RVs in the 
online tables. The new and old RVs are plotted in Figure QJ 
for WASP-22 and in Figure EJi for WASP-26, folded on the 
ephemerides of Table [3] 

We estimated the sky-projected stellar rotation velocity^ 
v sin / from the HARPS spectra by fitting the profiles of several 
unblended Fe i lines. For this we determined an instrumental 
broadening of 0.060 ± 0.005 A f rom telluric lines around 6300 A 
and used the Bruntt et al. l d2oToh calibration to assume values for 
macroturbulence: v mac = 3.2 + 0.3 km s _1 for WASP-22 and v mac 
= 3.0 ± 0.3 kms" 1 for WASP-26. We obtained v sin / values of 
4.5±0.4kms- 1 for WASP-22 and 3. 9 + 0.4 km s" 1 for WASP-26. 

We determined the typical activity indices (logR' HK and 
Smw) of each star by measuring the weak emission in the cores 
of the Ca n H+K lines in the HARPS spectra with the highest 
SNR per pixel at 550 nm: 13 spectra with SNR>3 2 for WASP- 
22 and 1 1 spectra with SNR>27 for WAS P-26 dNoves et al.l 
Il984t ISantos etal.ll2000t iBoisse et al.ll2009l) . The values given 
in TableQ~]are the averages and standard deviations of the values 
from individual spectra, determined using B - V values of 0.57 
for WASP-22 and 0.59 for WASP-26. For each star, we estimated 
the true stella r-rotation period, P mU using the activity-rotation 
calibration of Mamaiek & Hillenbrand (2008), and combined 
this with the stellar radius derived in Section [3] to estimate the 
true stellar-rotation velocity, v (Table [TJ. For each star, the sim- 
ilarity of the values of v and v sin / suggests the stellar spin axis 
is not significantly inclined relative to the sky plane. 



2.2. Photometric transits 

The cadence of spectroscopic transit observations tends to be 
much lower than that of photometric observations and, at least 
for HARPS measurements of WASP systems, the spectroscopic 
exposure time tends to be similar to the durations of transit 
ingress and egress. As such, the accurate measurement of A re- 
quires the times of the four transit contact points during the 
spectroscopic transit to be accurately determined. This is best 
achieved using high-quality photometric transit observations. As 
M10 presented only a partial transit of WASP-22b and as the 
two transits of WASP-26b presented by S10 were obtained nine 
months prior to the RM measurement, we obtained additional 
transit photometry of the two systems. 



2 We use / for the angle between the stellar spin axis and the sky 
plane and i for the angle between the planetary orbital axis and the sky- 
projection of the stellar spin axis. 
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Table 1. Stellar parameters 



Star 


T s g 




l°gg» 


[Fe/H] 


vsin/ 


log^ K 


5 MW 


^rot 


V 




(K) 




(cgs) 




(krns" 1 ) 






(days) 


(kms" 1 ) 


WASP-22 


6000 ± 


100 


4.5 ± 0.2 


0.05 ± 0.08" 


4.5 ±0.4 


-4.90 ± 0.06 


0.17 ±0.1 


14 ±3 


4.3 ±0.8 


WASP-26 


5950 ± 


100 


4.3 ± 0.2 


-0.02 ± 0.09 


3.9 ±0.4 


-4.98 ± 0.07 


0.16 + 0.1 


19 ±3 


3.5 ±0.5 



Notes. The values of T e $, logg, and [Fe/H] are from M10 for WASP-22 and from S10 for WASP-26. The other parameters are from this paper. 

(a) Due to a typographical error in M10, [Fe/H] was reported as -0.05 ± 0.08. 



2.2.1. WASP-22b 

We used TRAPPISTH a 60- cm robotic telescope located at ESO 
La Silla dGillon et aLll201 ll) to observe two transits of WASP- 
22b. Both transits were observed through a special I+z filter, in 
the 1x2 MHz read-out mode, and with lxl binning, resulting 
in a typical combined readout and overhead time of 8 s and a 
readout noise of 13.5 e~. With a pixel size of 15 x 15 //m, the 
pixel scale is 0.65" pixel -1 and the field of view is 22 x 22'. The 
guiding system of TRAPPIST, based on absolute astrometry per- 
formed on each image, kept the stars at the same positions on the 
chip to within 3" (5 pixels) over the course of the observations 
(Jehin et al.,2011). 

The data from both transits were processed in the same way. 
After a standard pre-reduction (bias, dark and flat-field correc- 
tion), the stellar fluxes were extracted from the i mages using th e 
IRAF/DAOPHOlQ aperture photometry software dStetsonl [l987). 
We tested several sets of reduction parameters and chose the set 
that gave the most precise photometry for the stars of similar 
brightness to WASP-22. We carefully selected a set of reference 
stars and then performed differential photometry. 

The first transit of WASP-22b that we observed occurred on 
2010 Nov 18. We obtained little pre-transit data as the transit 
started soon after dusk. We also obtained little post-transit data 
as WASP-22 passed through the meridian soon after the transit 
ended, at which time a meridian flip would have been necessary 
as TRAPPIST uses a German equatorial mount. The post-flip 
light curve would thus have been of little value as, due to po- 
tential systematics, it would have needed to be detrended sepa- 
rately to the post-flip light curve. The second WASP-22b transit 
we observed occurred on 2010 Dec 2. We needed to perform 
a meridian flip shortly after beginning observations and shortly 
before the transit began, so we divided the data into pre-flip and 
post-flip light curves. 

We also observed the 2010 Dec 2 transit of WASP-22b 
with EulerCam on the 1.2-m Swiss telescope at ESO La Silla. 
Installed at the Swiss telescope in September 2010, EulerCam 
uses a 4k x 4k, back-illuminated, deep-depletion e2V silicon 
chip. Read-out of the entire chip takes ~6 s in four-port mode 
and ~25 s in one-port mode. In both modes the readout noise is 5 
electrons. The field of view of EulerCam is 15.6 x 15.6', though 
the corners are shadowed by the filter wheel. With a pixel size of 
15 x 15 jjm the pixel scale is 0.23" pixeL 1 . EulerCam employs 
a very similar guiding system to that of TRAPPIST. Details of 
EulerCam will be provided in Lendl et al. (in preparation). 

We observed WASP-22 for 4.8 hr through a Gunn r filter 
with slight defocus and using one-port readout (as there were is- 
sues with four-port readout at the time). Conditions were good, 

3 TRAnsiting Planets and Planeteslmals Small Telescope; 
http://arachnos.astro.ulg.ac.be/Sci/Trappist 

4 IRAF is distributed by the National Optical Astronomy 
Observatory, which is operated by the Association of Universities 
for Research in Astronomy, Inc., under cooperative agreement with the 
National Science Foundation. 



with a seeing of 0.5-1.0". The target remained below airmass 
1.8 as we observed from 35 min prior to the start of transit un- 
til 45 min after the transit ended. We bias-subtracted and flat- 
fielded the images before performing differential aperture pho- 
tometry using a reference source created by co-adding the two 
most favourable reference stars. 

The two transit light curves of WASP-22b from TRAPPIST 
and the one from Euler are shown in Figure QJ. We also incor- 
porated in our analysis of WASP-22 the photometry of M10: 
two seasons of WASP photometry covering the full orbit and a 
partial transit measured by Faulkes Telescope South through a 
Pan-STARRS z filter. 

2.2.2. WASP-26b 

We observed a full transit of WASP-26b on 2010 Oct 28/29 
through a Sloan r' filter using the SBIG STL- 100 IE CCD 
camera installed at the prime focus of the 0.35-m Celestron 
C14 Schmidt-Cassegrain telescope at the robotic Oversky 
Observatory, La Palma. The field of view is 19.9 x 19.9' and 
the exposure time was 15 s, with only a small degree of defocus 
employed so as to avoid blending with a star located 15" from 
the target. The data were processed using MUNIPACK0 and the 
transit light curve is shown in Figure [2^. 

We also incorporated in our analysis of WASP-26 the pho- 
tometry of S10: two seasons of WASP photometry covering the 
full orbit and two full transits measured through a Pan-STARRS 
z filter, one by Faulkes Telescope North (FTN) and the other by 
Faulkes Telescope South (FTS). In-transit noise in the FTN light 
curve of S 10 (see their Figure 4) appears to have affected its nor- 
malisation. Thus, the transit depth and, in turn, the stellar and 
plantetary radii may have been biased. As such we chose to de- 
trend that light curve using only the post- transit data. S10 noted 
that WASP-26 is blended in the WASP data with the common- 
proper-motion star located 15" away that is ~2.5 mag fainter. 
We made an approximate correction to the WASP photometry 
for this contamination so as to prevent dilution of the transit. 

3. Combined model and analysis 

For each system, we determined the system parameters from a 
simultaneous fit to the data described in the previous section. The 
fit was performed using the current versio n of the Markov-Chain 
Monte Car lo (MCMC) code desc ribed bv lCollier Cameron et al] 
(12007b and lPollacco et all (120081) . 

The transit ligh t curve s were modelled using the formulation 
of Mandel & Agoll (120021) with the assumption that the planet is 
much smaller than the star. Limb-darkening was accounted for 
using a four-coefficient, nonlinear limb-darkening model, using 
coefficients appropr iate to the passbands from the tabulations of 
IOarel(l2000l 12004). The coefficients were interpolated once us- 
ing the values of log g* and [Fe/H] from M10 and S 10 (TableQ}, 

5 http : //c-munipack . sourcef orge . net 
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Fig.l. (a): The CORALIE (blue circles) and HARPS (green Fig. 2. (a): The CORALIE (blue circles) and HARPS (brown 

diamonds) RVs of WASP-22 from M10, and the HARPS RVs diamonds) RVs of WASP-26. (b), (c) and (d): Expansions of the 

from this paper (brown diamonds), (b): An expansion of the re- region around the spectroscopic transit, as measured by HARPS 

gion around the spectroscopic transit, as measured by HARPS on 2010 Sep 12, in which the RM effect is not modelled (b), 

on 2009 Nov 19. (c): The transit light curves, from top to bot- the RM effect is modelled and a prior is imposed on vsin/ (c), 

torn, obtained by TRAPPIST on 2010 Nov 18 and 2010 Dec 2 the RM effect is modelled without a prior on v sin / (d). (e): The 

through an I+z filter and by EulerCam on 2010 Dec 2 through Sloan r'-band transit light curve obtained at Oversky on 2010 

a Gunn r filter. The data obtained by TRAPPIST on 2010 Dec Oct 28/29. 
2 prior to the meridian flip are denoted by green triangles and 
blue dots denote the data obtained post-flip. In each panel, the 

data are phase-folded on the ephemeris given in Tableland the rameterised by the epoch of mid-transit 7q, the orbital period P, 

best-fitting models of Section0are superimposed. the planet-to-star area ratio (R p i/R„) 2 , the duration of the tran- 
sit from initial to final contact T14, and the impact parameter 
b = a cos (the distance, in fractional stellar radii, of the 

but were interpolated at each MCMC step using that step's value transit chord from the star's centre in the case of a circular or- 

of r e ff . The coefficient values corresponding to the best-fitting bit), where a is the semimajor axis and i is the inclination of the 

value of T e ff are given in Table [2] The transit light curve is pa- orbital plane with respect to the sky plane. At each MCMC step, 
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each light curve was decorrelated with a linear function of phase 
using singular value deconvolution. As noted in Section [2721 we 
detrended the FTN transit of WASP-26b from S10 using only 
the post-transit data, specifically those measurements with an or- 
bital phase greater than 1.018. From a fit to all data, this gives 
a slightly smaller {R p \/R«) 2 of 0.01022 ± 0.00034 as compared 
to the 0.01075 ± 0.00033 obtained when detrending using the 
entire light curve. This results in slightly smaller derived stellar 
and planetary radii (1.303 + 0.059 R Q c.f. 1.334 + 0.051 R Q and 
1.281 + 0.075 flj U p c.f. 1.346 ± 0.067 i?j U p). 

The eccentric Keplerian radial-velocity orbit is parame- 
terised by the stellar reflex velocity semi-amplitude K\ , the sys- 
temic velocity y, an instrumental offset between the HARPS 
and CORALIE spectrographs A^harps, a linear drift in the sys- 
temic velocity y, and y[e cos a> and y[e sin a> (Collier Cameron, 
in preparation), where e is orbital eccentricity and u> is the 
argument of perias tron. The RM e ffect was modelled using 
the formulation of iGimene 1 J2006h and is parameterised by 
Vv sin/ cos A and Vvsin/ sin A. 

The linear scale of the system depends on the orbital sepa- 
ration a which, through Kepler's third law, depends on the stel- 
lar mass M*. At each step in the Markov chain, the latest val- 
ues of p„, T e g and [Fe/Hl are in put in to the empirical mass 
calibration of Eno ch et alJ (1201 Oh to obtain M*. The shapes of 
the transit light curves and the ra dial-velocity curve constrain 
p„ dSeager & Mallen-Or nelas 2003), which combines with M* 
to give R*. T e ff and [Fe/H] are proposal parameters constrained 
by Gaussian priors, by means of Bayesian penalties on^ 2 , with 
mean values and variances derived directly from the stellar spec- 
tra (TableHJ. 

As the planet-star area ratio is determined from the measured 
transit depth, R p i follows from /?». The planet mass M p \ is calcu- 
lated from the values of K\ and M„ and planetary density p p \ and 
surface gravity log g p \ then follow. We also calculate the ingress 
and egress durations, T[2 and T34, and the planetary equilibrium 
temperature r eq i, assuming zero albedo and efficient redistribu- 
tion of heat from the planet's presumed permanent day side to 
its night side. 

At each step in the MCMC procedure, model transit light 
curves and radial velocity curves are computed from the pro- 
posal parameter values, which are perturbed from the previous 
values by a small, random amount. The y 2 statistic is used to 
judge the goodness of fit of these models to the data and a 
step is accepted if x 2 is lower than for the previous step. A 
step with higher x 2 is accepted with a probability proportional 
to exp(-Ay 2 /2), which gives the procedure some robustness 
against local minima and leads to the thorough exploration of 
the parameter space around the best-fitting solution. To give 
proper weighting to each photometry dataset, the uncertainties 
were scaled at the start of the MCMC so as to obtain a pho- 
tometric reduced-^ 2 of unity. With a similar purpose, a jitter 
term of 4 ms _1 was added in quadrature to the formal errors 
of the HARPS RVs of WASP-22, as may be due to stellar ac- 
tivity (though there is no evidence of this here), correlated noise 
or the finite number of data-points. This was determined from 
an initial MCMC from which the in-transit RVs were excluded. 
With spectroscopic reduced-^ 2 values of less than unity, it was 
not necessary to add any jitter to the RVs of WASP-26. 

To prevent a specific stellar activity level during a spectro- 
scopic transit observation from biasing the fitting of the RM ef- 
fect, RVs through a transit have previously been grouped sep - 
arately to RVs spread over a long baseline (Tri aud et al.l l2010). 
WASP-22 and WASP-26 both appear to be chromospherically 
inactive though (M10; S10), so we grouped RVs only by spec- 



trograph so as to be able to determine the rate of drift in the sys- 
temic velocity whilst allowing for an instrumental offset. We did 
check that this choice did not affect the spin-orbit determination. 
The systemic velocity at the time of the spectroscopic transit is 
set by the RVs taken around that time, thus assisting in the accu- 
rate measurement of the projected spin-orbit angle. In the case of 
WASP-22, combining these HARPS RVs with those measured a 
year prior results in an improved determination of the linear drift 
in the systemic velocity noted in M10 and in tighter constraints 
on orbital eccentricity. 

For each system, the improvement in the fit to the RV data 
resulting from the use of an eccentric orbit model is small and 
is consistent with the underlying orbit being circular. We thus 
impose circular orbits, as is prudent to do for such short-period, 
Jupit er-mass planets in the absence of evidence to the contrary 
fe.g. lAnderson et aTll201 ll) . This did not result in significant dif- 
ferences in the best-fitting values of the system parameters or 
their associated uncertainties. The 3-cr upper limits on e are 
0.063 for WASP-22b and 0.050 for WASP-26b. 

We imposed a Gaussian prior on vsin/, by means of a 
Bayesian penalty on^ 2 , with mean and variance as determined 
from the HARPS spectra (see Section |2J]>. For WASP-22, this 
prior prevents v sin / from wandering to unrealistically large val- 
ues of up to 168 kms -1 . This occurs due to a degeneracy be- 
tween v sin / and A when b is low: an RM effect with a similar 
amplitude and shape can result both when a planet transits near- 
parallel to the equator of a slowly rotating star and when a planet 
in a near-polar orbit transits a rapidly rotating star. The MCMC 
posterior distribution of vsin/ and A shown in Figure [3^ results 
when a prior is imposed, and the distribution shown in Figure[3j3 
results when no prior is imposed. 

We checked if the fitting of the RM effect, which involves 
the inclusion of the two parameters vsin / and A, results in a 
significant improvement to the fit to the data for each system. 
For this we performed an F-test of the null hypothesis that we 
did not detect the RM effect, using a probability of P(F) = 0.05 
as our significance threshold. 

For WASP-22, we find a highly significant improvement to 
the fit when modelling the RM effect with a prior on vsin/: 
P(F) ~ 0. So, we have made a significant detection of the RM ef- 
fect of the WASP-22 system, with A = 22 ± 16°, which is clearly 
visible in the star's RVs (FigureQJi). Without the quadrature ad- 
dition of the 4 ms _1 jitter term to the HARPS RVs, a value of 
A-25 + 13° is obtained. 

For WASP-26, we are unable to reject the null hypothesis as 
the data are fit equally well by both a model with no RM ef- 
fect and by a model with the RM effect and a prior on v sin /, 
for which A = 6 ± 6° (P(F) ~ 1; see Figures QJ> and [[};). 
When modelling the RM effect with no prior on v sin /, for which 
A = 7 ± 50°, we find the improvement in the fit is too slight 
to reject the null hypothesis (P(F) = 0.50; Figure QJl). Also, 
the resulting best-fitting projected stellar rotation speed (vsin/ 
= 0-7^H kms -1 ) is discrepant with the spectral determination 
(vsin/ = 3.9 ± 0.4 kms" 1 ). It can be seen in Figure|2]; that two 
in-transit RVs lie above the model. The seeing was poorer when 
those two measurements were obtained (1.1") than during the 
rest of the RM sequence (0.5-0.8"). We explored excluding the 
four RVs (the two in-transit and two others outside of transit) 
with seeing greater than 1", but found that a model with the RM 
effect and a prior on vsin/, for which A = 13 ± 6°, was still not 
strongly favoured over a model with no RM effect: P(F) = 0.09. 
On the night of the RM measurement, three RVs were obtained 
prior to the re-focussing of the instrument. These data-points lie 
below the best-fitting model (Figure 0. Having omitted these 
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Fig. 3. The MCMC posterior distributions of A and v sin / for 
WASP-22 when (a) a Gaussian prior of 4.5 + 0.4 kms~'is im- 
posed on v sin / and (b) no prior is imposed on v sin /. Note that, 
though the abscissa ranges of each graph are equal, the ordinate 
ranges are not. The white, blue and black contours are, respec- 
tively, the 1-, 2- and 3-cr confidence limits. The shading of each 
bin is proportional to the logarithm of the number of MCMC 
steps within. Aligned spin-orbit axes would have A = 0°, which 
is marked by the dotted, brown line. 



three RVs, a similar (slightly worse) fit to the RVs data resulted 
when fitting the RM effect with the prior on v sin / as compared 
to when not fitting the RM effect. As such, we did not make 
a significant detection of the RM effect for WASP-26 and we 
adopt a Keplerian orbit with no RM effect as our model for that 
system. We note, though, that in Figure \2)p there appears to be 
some structure in the in-transit RVs. 

We explored the degree to which the new transit photometry 
was assisting in the fitting of the RM effect by excluding it from 
MCMCs for the two systems. For WASP-22b, the uncertainty in 
the transit mid-point at the epoch of the RM measurement was 
reduced from 3 to 1 minutes, and the uncertainty in the durations 
of the transit, ingress and egress fell from 6.5 to 1.4 minutes. 
For comparison, the ingress and egress durations are 20 minutes. 
Also reduced were the uncertainties in b (from 0. 16 to 0. 12) and 
Rpi/R* (from 0.024 to 0.013). These parameters, relevant to the 
fitting of the RM effect, and others were previously less certain 
due to the lack of available high-quality transit photometry. This 



led to the imposition of a main-sequence prior on the host star 
in M10, which was not imposed here in the MCMC omitting 
the new photometry. Without the new photometry A = 6^° i s 
obtained. For WASP-26b, the uncertainty in the transit mid-point 
at the epoch of the RM measurement was reduced from 2 to 
1 minutes, and the uncertainty in the durations of the transit, 
ingress and egress fell from 7.7 to 1.4 minutes. The ingress and 
egress durations are 34 minutes. 

By measuring the inclination of the stellar spin axis with re- 
spect to the sky plane /, we could determine the true, rather than 
the sky-projected, obliquity. A measurement of the stellar rota- 
tion period P mt , in combination with R*, would give the stellar 
rotation velocity v. By comparing this to the sky-projected ro- 
tation velocity vsin/, determined from the spectra, we could 
obtain /. P rot can be measured if a star has an active photo- 
sphere that induce s rotational modulation in its light curve (e.g. 
Helli er et al 1 l201lh . However, M10 found no evidence for mod- 
ulation in the WASP-22 light curves, and S10 found none in the 
WASP-26 light curves. 

3.1. Results and discussion 

The results of the MCMC model fits are presented in Table [3] 
The corresponding transit and orbit models are superimposed 
on the new transit light curves and on all RVs, respectively, in 
Figure[TJfor WASP-22 and in Figure|2]for WASP-26. 

With A -22 + 16°, WASP-22b is in a prograde orbit around 
its host star. The best-fitting model suggests that the planetary or- 
bital axis is slightly misaligned with the stellar spin axis, but the 
current uncertainty allows for quite a range in A (Figure^). The 
similarity of our estimat es of v sin / and v (T able[T]i, which agrees 
with the prediction of Schlaufma^ (1201 Oh that v - 3.2 ± 1.0 
kms -1 , suggests that the spin axis of WASP-22 is not signif- 
icantly inclined relative to the sky plane (i.e. ^ ~ A). With 
T eS = 6000 ± 100 K, WA SP-22 is consistent with the observa- 
tion of IWinn et al.l d2010al) that the orbits of planets around stars 
cooler than ~6250 K tend to be aligned. Conversely, the orbits 
of pl anets around hotter stars tend to be m isaligned, as noted by 
both lWinn et all d2010al) and ISchlaufm an (2010). A few possi- 
ble expections to this trend are now known, including HAT-P- 9 
with r eff = 6350 K and A = -16 ± 8° dMoutou et : alJl201 ll) . 
Though it may be t hat such systems were never very misaligned. 
Win n et al.l d2010al) suggested the observed trend either indicates 
that planetary formation and migration mechanisms depend on 
stellar mass, or that the larger convective envelopes of cooler 
stars are able to more effectively realign with orbits via tidal 
dissipation, though this does depend on speculative, long-lived 
core-envelope decoupling. 

M10 detected a linear trend of y — 40 + 5 ms _1 in the RVs 
of WASP-22, which span a duration of 16 months. The new 
HARPS RVs presented herein make possible a more precise de- 
termination of the trend of j = 37.8 ± 2.2 ms~' (Figure 0). 
We can obtain an order of magnitude estimate of the mini- 
mum mass of the third body, M3 sin £3, scaled by its orbita l dis- 
tance, 03, from the primary. Following IWinn et al.1 ((2010b), by 
assuming the third body to be in a near-circular orbit and to 
have a mass much smaller than that of the primary, we may set 
y ~ GMi, sin 13 /o^ to obtain: 



M3 sin h 



M 



J up 



\5AU/ 



5.3 ±0.3. 



(1) 



The RM effect of WASP-26b, predicted to have a semi- 
amplitude of ~15 ms~', was not detected with significance in 
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our data, for which the median in-transit uncertainty was 7ms. 
This was due to observation noise and the small amplitude of the 
signal. For comparison, we measured the semi-amplitude of the 
RM effect of WASP-22b to be -30 ms" 1 , which agrees well 
with our prediction of 27 ms _1 . The median in-transit RV un- 
certainty was 8 ms~'. The predicted amplitude of the RM ef- 
fect of WASP-26b is smaller than that of WASP-22b due to the 
higher impact parameter and, to a lesser extent, the slower stel- 
lar rotation velocity. The similarity of our estimate s of vsin/ 
and v ( Tabled}, which agrees with the prediction of Schlaufman 
(2010) that v — 3.1 ± 0.3 kms~', suggests that the spin axis of 
WASP-26 is not significantly inclined relative to the sky plane. 
Thus it is unlikely that the non-detection of the RM effect is due 
to the orientation of WASP-26 being close to pole-on. Though it 
is possible that the planet is in a polar orbit, it can not be transit- 
ing over the stellar spin axis as the impact parameter is non-zero. 
Hence, there will be an observable RM signal and thus further 
measurements are warranted. 

The available RVs of WASP-26 are not a very sensitive probe 
of a long-term drift as the CORALIE RVs span only 95 days 
and the HARPS RVs span only 5 days, with a gap of 1 year 
between the two datasets. Any drift could therefore have been 
taken up within the fitted instrumental offset between the two 
spectrographs. However, the similarity in the offsets for the two 
systems (Tabled suggests that any drift in the systemic velocity 
of WASP-26 is likely to be smaller than that of WASP-22. 
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Fig. 4. Radial velocities of WASP-22 showing a linear drift in 
the systemic velocity of y — 37.8 ± 2.2 ms _I yr _1 . The colour 
and symbol key is the same as in Figure [1] The Keplerian orbit 
of Table[3]was subtracted. The value of y from Table[3]is plotted 
as a solid line, relative to JD = 2 455 013, which is the centre-of- 
mass of the RV data, weighted by the square of the measurement 
precision. The drift's 1-cr error bars are plotted as dashed lines. 
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Table 2. Limb-darkening coefficients 



System 


Instrument 


Observation band 


Claret band 


a { 


a 2 


«3 


«4 


WASP-22 


WASP/Euler 


Broad (400-700 nm) / Gunn r 


Cousins R 


0.578 


-0.049 


0.494 


-0.294 


WASP-22 


TRAPPIST 


Cousins /+Sloan z' 


Sloan z' 


0.652 


-0.348 


0.638 


-0.326 


WASP-26 


WASP 


Broad (400-700 nm) 


Cousins R 


0.573 


-0.027 


0.464 


-0.281 


WASP-26 


FTN / FTS 


Pan-STARRS z 


Sloan z' 


0.651 


-0.337 


0.624 


-0.321 


WASP-26 


Oversky 


Sloan r' 


Sloan r' 


0.551 


0.048 


0.415 


-0.266 



Table 3. System parameters from RV and transit data from our adopted, circular solutions 



Parameter 



Symbol 



WASP-22 



WASP-26 



Unit 



Orbital period 

Epoch of mid-transit (HJD,UTC) 
Transit duration 

Transit ingress (= egress) duration 
Planet-star area ratio 
Impact parameter 
Orbital inclination 



P 

T 
Tu 

Tl2 = 7^34 

(tfp/R,) 2 
b 
i 



Stellar reflex velocity semi-amplitude K { 

Semi-major axis a 

Offset between HARPS and CORALIE Ay H ARPS 

Systemic velocity at time T y 

Linear drift in systemic velocity y 

Orbital eccentricity e 

Sky-projected spin-orbit angle A 



Sky-projected stellar rotation velocity 


vsin/ 




v sin / cos A 




v sin I sin A 


Stellar mass 


M, 


Stellar radius 


R, 


Stellar surface gravity 


logg, 


Stellar density 


p, 


Stellar effective temperature 




Stellar metallicity 


[Fe/H] 


Planetary mass 


M pl 


Planetary radius 


^ P i 


Planetary surface gravity 


l°ggpl 


Planetary density 


P P i 


Planetary equilibrium temperature 





3.5327313 ± 0.0000058 
2455497.39967 ± 0.00025 
14501+ - 00114 
01369 + ™ 24 

OJ]0074 

0.00954 ± 0.00018 

0. 25 ±0.12 
88.26 ±0.91 

73.2 ± 1.7 
0.04698 ± 0.00037 
17.60 ± 0.80 
-7 187.3 ±3.5 
37.8 ±2.2 

(adopted) 
< 0.063 (3 cr) 

22 ± 16 
4.42 ± 0.34 
3.98 ± 0.43 
1.63 ± 1.13 

1.109 ±0.026 

1 tiq+0.052 

1. z.l5>_ 0033 

4 t 1 n+0.020 

612+™ 
5958 ± 98 
0.050 ± 0.080 



0.588 
1 



: 0.017 



| cq+0.061 

3.001^1 
0.378!»f 
1466 ± 34 



2.7566004 ± 0.0000067 d 
2455228.38842 ± 0.00058 d 
0.0993 ± 0.0018 d 
0.0239 ± 0.0026 d 
0.01022 ± 0.00034 
0.812 ±0.020 
82.91 ± 0.46 



137.7 ± 1.5 
0.03985 ± 0.00033 

20.62 ±0.14 
8 459.297 ±0.073 

(adopted) 
< 0.048 (3 cr) 



1.111 ±0.028 
1.303 ±0.059 
4.253 ± 0.034 
0.502 ± 0.062 

5939 ± 100 
-0.020 ± 0.091 

1.028 ± 0.021 
1.281 ± 0.075 
3.156 ±0.048 
0.488 ± 0.082 
1637 ±45 



ms 
AU 

ms" 1 
ms" 1 
ms^'yr -1 



kms~' 
kms -1 
kms -1 

M 
Re 
(cgs) 

Po 
K 



M Su] , 

(cgs) 

pjup 
K 
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